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Q\ ; ABSTRACT 

0^ ' We developed theory of particle acceleration inside the relativistic rotating electron- 

positron force-free jet with spiral magnetic field. We considered perturbation of sta- 
^ ' tionary magnetic field structure and found that acceleration takes place in the regions 

where the Alfven resonant condition with the eigenmodes in the jet is fulfilled, i.e. 
^< ' where the local Alfven speed is equal to the phase speed of an eigenmode. Accelera- 

tion process and synchrotron losses combined together form power law energy spec- 
trum of ultrarelativistic electrons and positrons with index between 2 and 3 depending 
^ ^ upon initial energy of injected particles. Synchrotron emission of these electrons and 

\^ • positrons in spiral magnetic field of rotating force-free jet has been calculated. Polar- 

T^lj- I ization properties of the radiation has been obtained and compared with existing VLBI 

' polarization measurements of parsec-scale jets in BL Lac sources and quasars. Our 

I results give natural explanation of observed bimodality in alignment between electric 

field vector and jet axis. Degree of polarization and velocity of observed proper motion 
of bright knots depend upon angular rotational velocity of the jet. Thus, comparing 
them to each other, we can estimate angular rotational velocity in jets. These results 
indicate that the fact, that generally in BL Lac objects electric field vector is oriented 
parallel to the jet axis while in quasars perpendicular to the jet axis, may be due to 
O I intrinsically larger angular rotational velocity and large winding of magnetic field in 

^ ■ BL Lac jets than in quasar jets. 

Key words: radiation mechanisms: nonthermal - magnetic fields - plasmas - galax- 
ies :jets. 
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1 INTRODUCTION 

The nature of the Active Galactic Nuclei (AGN) is still unclear. The most common point of view is that there is a supermassive 
black hole in the centre of active galaxy with mass approximately 10* W^Mq (Rees, 1984). Accretion on the black hole leads 
to the formation of powerful turbulent processes in accretion discs which effectively heat the plasma and generate magnetic 
field of the order of 10*G. Ejection of plasma from discs with frozen magnetic field and emission lead to the formation of 
colhmated streams which transfer energy to large distances, up to 10^ pc. Besides, the processes of generation of electron- 
positron plasma may occur in the magnetospheres of the black holes as well as in the vicinity of accretion discs. Particles 
are created in collisions of high-energy photons produced by inverse Compton scattering of ultraviolet photons emitted from 
disk with fast particles accelerated by magnetic reconnection within disk or within magnetospheres of the rotating black 
hole, where equilibrium charge density changes its sign (see, for example Beskin, Istomin & Pariev, 1992b). Electron-positron 
plasma forms pinched streams (jets) which have radii of the order of one parsec. The energy density of electromagnetic field 
in these jets greatly exceeds the energy density of electrons and positrons taking into account their rest energy. Rotation of 
the black hole and the accretion disc is transmitted to the jet due to magnetic coupling, where strong radial electric field Er 
arises. Drift in this field and in the longitudinal magnetic field Bz frozen into plasma leads to differential rotation. Electric 
current transferred by the jet produces azimuthal magnetic field -B,^, and magnetic field lines are twisted into spirals. Energy 
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in the jet is transferred by the Pointing vector which is proportional to the product of radial electric field and azimuthal 
magnetic field. Therefore, there is no problem of energy transfer along the jet on large distances as there is in the case when 
the main energy is in the hydrodynamic motion. 

The most amazing facts of the jets from AGN are their high coUimation (the ratio of the length to the diameter is 
a few orders of magnitude), stability and apparent superluminal velocities of distinct knots (Begelman, Blandford & Rees, 
1984). The hydrodynamic and MHD stability of jets was investigated in many works (Blandford & Pringle, 1976; Torricelli- 
Ciamponi & Petrini, 1990; Appl, Camenzind, 1992). As to the stability of the electron-positron jets, rotating and moving with 
relativistic speed, it was shown in papers by Istomin & Pariev, 1994 and Istomin & Pariev, 1996 that they are stable with 
respect to axisymmetric as well as spiral perturbations. The physical reason for stability is the shear of the magnetic field, 
which in the case of small hydrodynamic pressure of plasma stabilizes small oscillations. Perturbations do not increase with 
time (Imcj = 0) or have small damping decrement {Imu} « 10~^Recij) because of the resonance with Alfven waves lo' = fcyC. 
In electron-positron plasma, when the specific energy density and pressure are much less than the energy density of magnetic 
field, Alfven velocity is equal to the speed of light in vacuum c (w' and k' are the frequency and the wave vector in the 
plasma rest frame). Resonance condition is fulfilled on the specific magnetic surface and, in the case of cylindrical jet, at the 
definite distance from the axis. In the vicinity of that surface magnetic and electric fields of the wave reach large magnitudes. 
The particles of the jet are accelerated by the electric field here, which leads to the absorption of the energy of perturbation. 
Thus, the stability of the jet is connected directly with the production of energetic particles in the current. This resolves the 
problem of in situ acceleration of electrons and positrons producing synchrotron emission from knots observed along jets. In 
fact, as it is known (see Begelman at al., 1984), energetic particles accelerated in central region can not penetrate in the jet, 
they must loose their transverse momentum by synchrotron radiation already in the basement of the jet. Characteristic time 
of synchrotron losses is given by 

_ 3 C7|| 
4 rew^o 

Here 7|| is the longitudinal Lorentz-factor of the particle, 7|| = (1 — ii^/c^)"^/^; re is the classical radius of the electron 
Te = /mc^; ljbo = eB/mc is the electron cyclotron frequency. Given 7|| = 10, B = lO^G radiation time Tg is approximately 
3 • 10® sec whereas the typical time of the flyby through the region of strong fields r/ w i/c w lO^sec. Here I is the diameter 
of the jet, £ w 3 ■ lO^'^cm. 

In this paper we assume that the energy transferred along the jet, is mainly the electromagnetic energy which can be 
transmitted for long distances along the jets as along the wires. In this case it is very possible that part of that energy is in 
waves propagating along the jet, which are eigenmodes of cylindrical beam. The source of the wave motion is nonstationary 
processes in the magnetospheres of black hole and accretion disc. Short time variability on scales from days to months 
is actually observed in AGN (Mushotzky, Done & Pound, 1993; Witzel et al., 1993). Such a picture also gives a natural 
explanation of the superluminal motions but requires neither special orientation of the jet nor relativistic jet with high speed 
(/? = u/c > V2/2). As it was shown in papers by Istomin & Pariev, 1994 and Istomin & Pariev, 1996 there exist so called 
standing modes among eigenmodes in jets {vgroup = 0) which are not propagate with finite velocity along jet but are only 
subject to diffuse spreading and, therefore, their amplitude is maximal. Phase velocity of these modes is greater than the speed 
of light. Crests of the wave move along the jet with superluminal velocity causing acceleration of the particles on the Alfven 
surface. These regions can be the objects which are seen as bright knots with typical sizes of the order of the wavelength of 
the standing wave which is about the radius of the jet. 

The purpose of our paper is the consideration of the acceleration processes near the Alfven resonant surface, generation 
of the spectra of accelerated particles and calculation of the polarization of synchrotron radiation generated by these particles. 



2 ACCELERATION OF PARTICLES 

The equations of motion near the Alfven surface are rather complicated. Nevertheless, we can use drift equations (Sivukhin, 
1965) because the Larmor radius of relativistic electrons and positrons rc = S/eB « lO^cm (if we take B w 10"^ and 7 w 10) 
is much less then radius of the jet and the width of resonant surface ro (see expression (31*)). 

^ = -II I + X B] + i|[B X (BV)B] + Ji[B X VB]; 

dE _ dr Ev\ dB , , 

dt 

Here the electric field is assumed to be small compared with the magnetic one, £ is the energy of the particle, Jx = p\/B. 
Electromagnetic fields are equal to the sum of stationary fields (without subscript) and wave fields (subscript 1). For cylindrical 
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jet the stationary configuration of fields is (Istomin & Pariev, 1994) (c=l): 

B = B^e^ + B^e^- B^ = n'^rB^; 

u = KB + n^re^; (2) 



E = -fi^rfe^xB], 



where r, 0, z are cylindrical coordinates. In components 

B = (0,n^r,l)S,, 

E = (-n^r,0,0)S,. 

Here and K are functions of r, B^ does not depend on r. Stationary electric field, which has an absolute value of Q.^rBz, 
is not small compared to the magnetic field when OFr > 1, therefore we must consider our equations in the frame moving 
with plasma where E = 0. There are a lot of such reference frames because there is an arbitrary parameter K{r) in the 
expression for the velocity of plasma which determines the radial profile of longitudinal velocity. K(r)B^ is the velocity along 
the magnetic field which is not related to the rotation. We choose, however, the velocity u which minimizes the kinetic energy 
of the plasma in the stationary reference frame. 

_ (o,n^r,-(n^r)^) 

It seems that the plasma moves with that velocity in real jets. Fields in the wave calculated in terms of Lagrangian displacement 
^ are equal (here we drop out the phase coefiicient exp(— iwt + ikz + irrKf))) 

Bri = iB^Fi; 
B^, = -B.in'-rf + + ^T?); 

E., = B.in-r§.+^lin-r)-^Dy, 

-E^i = -iB^{uj - mfl^)^; 
E^i = iB^Q,^r{uj + k)i; 



where 



F = k + mQ^, 



S = uj — k — m / 



r 



For Lagrangian displacement ^ and dimensionless pressure disturbance p* = AirPi /B1 we have a system of differential equations 
obtained by Istomin & Paxiev, 1996: 



Here 



Ci = l {mn'" - (a^r)^(a; + k)) ; 
02 = —- 



U! + k 

C3 = -(a; + /c)(yl'-4n^'); 

A = k-uj + 2ma^ - (n^r)^(w + k). 
Let us expand equations (5) near the Alfven surface, i.e. near the point rA where A{rA) = 0. Denoting a; = r — we get 
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where A' = ^ 
ar 



. The general solution of these equations is 



a)=".(cr-%)-=(c:) 

where ai and 02 are constants. They are not arbitrary, rather they are fixed by the boundary conditions on ^ and p, obtained 
from solving general equation (5). We see that ^ and p* have a logarithmic singularity at the Alfven point. Later we adhere 
to our expansion over r — r a , assuming r = va, A = whenever the quantities under consideration have no singularities at 
r = rA- Let us calculate dr/dt and d£/dt in the first order by fields of the wave. Remembering that E = in our reference 
frame and knowing that terms with the energy £ contain small parameter which is the ratio of the Larmour radius to the jet 
radius or to the characteristic length of variations of zero order fields we can simplify drift equations as follows, 

' dr\ 

:W||e, 



\dt 



fdr\ Bi B(BjBJ [EixB] 



Here subscript denotes the zero order term and 1 denotes the first order term. The force of inertia is also present in the 

expression of dr/dt (sec Sivukhin, 1965), but the term with inertia force contains the mass, has the same order as terms with 
energy and is omitted consequently. For d£/dt we will also not take into account the inertia force because it is smaller than 
gradient terms such as fwj|[B x (BV)B]/ecB*. When we substitute Bi and Ei into (8) we get 

where 7 = -^/l + (fi^r)^ is the Lorentz factor corresponding to the velocity of the moving frame. ^ depends on z, 4>, t by 
exp{—iu!t+ikz+im<l)}. The phase of the wave ^ which is "seen" by the particle can be calculated using Lorentz transformation 
formulae >]/ — ~^t' + ^z' . Assuming dz/dt = we get >]/ — (t)|| — As we see the coefficient for the phase coincides with 

the coefficient in (8*). It confirms our calculations because ^ is the Lagrangian displacement. The question of the trajectory 
of the particle will be considered later. 

Now we proceed with calculating d£/dt. 

d£ ^ fdY\ £v±dBi 
After substitution of expressions for fields into (9) we get 

f = ^-i(-^)^(- + '=)[^-l. (10) 



or, combining with (8*) 

d£ _ £v\ f {u} + k)[F-, 



dr 1 — v\\ \ rS 



(11) 



As we see, £ is proportional to the r if we do not take into account dependence of the right hand side on £ Therefore it makes 
sense to consider motion of the particle "smoothed" over many periods. Taking the r-component of (8*), we have 

^ = ia;*$(x)e'-**, a;* = (t;,, - 1)^. (12) 

Near the Alfven resonance ^{x) has a logarithmic singularity 

dx 

— = u}*[A In X sm{ij* t) + B cos(u)*t)]. (13) 

It is clear that u)* A and liJ*B are dimensionless amplitudes of the perturbation (remember, that c=l). Let us average equa- 
tion (13) expanding it in small amplitude of oscillation of a particle in the field of the wave. We denote perturbations of the 
zero, first, second orders as xo, xi, X2- 

= u!* [A In Xo 8m{oj*t) + B cos(ai*t)] 
xi = — j41na;o cos(a;*i) + Bsin(a;*t) 
= a;*[Aln \xo — Alnxo cos(ui*t) + Bsm{ui*t) \ sin(a;*t) + B cos{uj*t)] 

dX2 _ 

dt 2xo 
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Now let us replace xo and X2 by x. We obtain smoothed equation of motion 



which has the solution 



dx _ « AB 



X = VuFABt. 



(14) 



Prom equations (11), (14) it follows that the particle gains energy drifting along the r axis. However , its acceleration 
rate is proportional to the transverse momentum which remains constant in accordance with the drift equations (1) 



dp± dJ± 

~dr °^ ~dr 



: 0. 



Here we do not take into account the change of the magnetic field connected with the wave field because of Bi <§; B. However, 

2 

because of synchrotron losses J± and p± are monotonically decreased according to (ws = u>bo — ^ — ) 



dJ± 

dt 



3 c 



(15) 



this leads to the situation whenpy > p±. Such anisotropic distribution can become unstable to the excitation of electromagnetic 
waves which results in isotropization of distribution function in the momentum space. We will treat instability of plasma 
according to Mikliaylovskii, 1968. Mikhaylovskii, 1968 showed that in the case > p± only perturbations with = must 
be taken into account. The frequencies of these perturbations with k _L B are the solutions of dispersion relation £33 — N'^ = 0, 
where £33 is the component of dielectric tensor along the magnetic field, N = ck/uj. The contribution of i sort of plasma into 
£33 is given by the general formula 



^33 



4ne 



\n= — oo 



oj — nujB 



dF 
dpz 



nuiB 



dF 
dpz 



Vz dF 
v± dp± 



(16) 



Here (. . .) means J . . . p±dp±dpz and the distribution function F is chosen so that J Fp±dp±dpz — no, no is the equilibrium 
particle density. From formula (16) it is seen that electron and positrons make an equal contribution to the £33 because the 
sum is symmetric to the simultaneous replacement of n with — n and of e with — e. After some calculations we get 



£33 - 1 = 



opz \ \ ujb J J v± dp± / \ ^-^ (jj^ — n^ujg Xujb J 



(17) 



Dispersion relation of the small oscillations has the form 



dpz ^ °' v± dp± 



(18) 



Study of (18) shows (Mikhaylovskii, 1968) that instabihty arises at small values of uj. Therefore let us assume a; to be small 
and find it 



C K 



dF 
dpj_ 



1 



(19) 



For further calculations we assume that kv± is of the order of ujb and estimate expression (19) taking into account only first 
terms of expansion of Bessel function over its arguments. 

, 2 



2 

U! = 



(20) 



Let us take integrals in (20) in parts remembering that ujb oc I/7. Here 7 = £ /mocp . As was shown in (Mikhaylovskii, 1968) 
the maximum increment is achieved at the values of ~ (oi^/vx). We will see later that instability at this n satisfies our 
requirements 



2 



7^ mi 



(21) 



Here Up = Anne^ /m is the plasma frequency squared. We continue our consideration by simplifying our assumption about the 
distribution function. Let all particles have equal p± and py . This is rather rough approximation but it allows us to understand 
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where the instabihty threshold is. So, remove all signs <> and for the sake of simplicity let mo = c = 1. 

^2 _ a|o 7V/3 + pj +p1- {pl +pt + 7'p^)/2 , . 

where /3 = 4immc^^/B'^ is the ratio of the particle pressure to the magnetic field pressure. Although /3 depends on energy let 
it to be constant to know where the threshold of the rcsonauco is in the case of fixed energy. If we consider ultrarelativistic 
particles p^ -C p**, we will have for the ratio of Pz to p± at which the imaginary part of u) arises 



pl _ 4-/3 + v^/3(17/3-8) 
Pl- 4(/3-l) • 

We denote the threshold of stability Pz/p± = a, o; > 1 (Expression (23) is not quite valid for very large /3 > 8 but we do 
not consider large (3 here). Im w arises when pz > ap±. As we see from (22) instability is fast because the increment of it is 
proportional to the cyclotron frequency. In our case it is the least typical time and it is less then the wave period by twelve 
orders of magnitude. Because of this fact we will always assume the distribution function to be such that Imco = 0. More 
concrete, we will consider Pz/p± = a = const. Assuming constant pitch angle for all particles may be confusing because as a 
result of the synchrotron emission of these particles we will obtain sharp maximum of the emission in directions close to this 
angle. But one has to remember that the emission must be integrated in the different parts of the jet, and always there will 
be particles found having velocities pointed to a given point in the sky, since in different parts of the jet we transit into a 
different reference frames. Thus the final result may be less rough than the intermediate. Now we understand that dJ±/dt ^ 
because of the anisotropic plasma instability and instead of the equation dJ± = constant we use Pz/p± = a = constant. For 
the derivative of energy we write equation (11) adding the synchrotron losses. 



d£ _ £vi iu; + k)[F-u>] 2 pi 



dt 1 — v\\ rS 3 mo 



Vx-- Wflofe- (24) 



Substituting Vx from (14) we obtain 



d£ _ £ {LO + k)[F -uj] f lo^AbY^^ 2 £' 2 



dt 1 + a2 - qVI + a2 rS V / 3 I + mo 

Let us introduce some conventional signs 

1 {uj + k)[F- uj] 



uJbo • (25) 



Q = 



1 -\- Q? — aVT+o^ S 
mo / {w + k)[F-c 



4 1 + 2a2 - 2aVl+a2 \ S J wl^rer'^'^* ' '■^^^ 

Here (w + k){F — ijj\/S is the fixed number for a given values of wavenumbers of the perturbation m and k. It is large near 
the fast magnctosonic resonance surface 5 = 0. The cocfRcient of \/{Lo%Qrer^uj*) is a large number . In the case of an AGN 
jet B « 10~^G (Begelman at al., 1984), r « Ipc it equals approximately to the 10'*. The value of w*^AB as we noted after 
formula (13) is the dimensionless amplitude of perturbation squared. Now taking into account our notations the solution for 
£ has the form 

Here £a = £{x = Q). Let us denote Qx/r = x' . We assume that x' can be large because of the quantity Q is large. It is 
seen that first £ exponentially increases and then decreases as 1/x' . We consider an asymptotic behaviour of £ in the case of 
different x and £o. In the case of £0 < £1 

-'^1 > ^^/^° (27*) 
x'e^ <£i/£o 



r £i/x', 

1 foe^'. 



In the case of £0 > £1 

c _ \ £i/x' , x >1 

~ \ 2fl/a;'^ x' < 1 

We reject the case when £0 > f 1 assuming that the typical initial energy is much less then £\ . 



3 FORMATION OF THE SPECTRUM OF ACCELERATED PARTICLES 

Our goal is to calculate the distribution function averaged over x' knowing the distribution function at a; = 0. We will consider 
the stationary case, when df /dt = 0. Let the distribution function to be given at the point x'o as a function of £' = £ /£i. 
Then the number of particles 

dN = f{£'o,x'o)d£'odx'o = f {£'o{£', x),x'o{£', x')) ^§r^d£'dx'. 
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Further we omit primes. Thus if we know the trajectory £o{£, x),xo{£, x) we may obtain the distribution function at the point 
£,x. In our case the Jacobian D{£q,x'o)/D{£',x') is not equal to unity because variables £,x are not canonical. So, 

F{£,x) = f(£oi£,x),xo{£,x))^^. (29) 

The formula analogous to the (27) but with the arbitrary reference point xo is the following 

£ = {x-l + {l/£{xo) -xo + 1) exp(-a; + xo)}~^ (30) 
Substituting £o expressed from (30) into (29) we obtain 

Fie.., . / ({-1 + .0 + (i - . + l) exp,. - .,)}-,..) (§^ - . ,31, 

/ ^|-1 + XO + - X- + 1^ cxp(x - X())| ,xo^ ^ exp(x - xq)- 



X 
X() 



We know that initial fast particles are born near the point a: = , therefore in the last formula we let xo = everywhere 
except Xo in the denominator. Physically f{£o,0) is the particle density at zero point and xo is the distance at which our 
formulae, for instance (14), do not work. Because of the small imaginary part of cj which is an attenuation increment of the 
wave (Istomin & Pariev, 1996) the divergence Inx at the resonance is cut at the point 

A{r = ta) 



ro « ,r;'~ , ■ (31*) 
dA/dr{r = rA) 

Here A{r = va) is not equal to zero as we assumed in the previous text, but determined by the small value of the imaginary 
part of uj. Thus, xo is the width of the resonance xo = Qro/rA- Eventually, 



F(f,x) = /o(^{-l+Q-x + l)exp(x)} f|expW^ 



(32) 

Xo 

Because the width of the resonance ro and the width of the acceleration region Ar , which is some function of ro , are small 
compared with the jet radius we average the distribution function over the jet radius near the resonance surface 

oo 

fi£) = ^ I Fi£,x) dx, Ax = ^Ar. (34) 



Ax 7 ' ' I'A 

Denoting N(£) = f(£)Ax we obtain 



oo 

Ni£) = J /o(^{-l+(i-x + l)exp(x)| 



-1\ g2 ^ 

-;|exp(x) — dx. (35) 

O Xo 





Let us change variable of integration to £o. The dependence of £o versus x under fixed £ is presented in Fig. 1. 
The value of £min is 

£rmn = l/(-l + eXp(l/£:)). (36*) 

Because the function x = x{£o) has two branches the equation (35) is transformed into 

s 

=xo J ^^^'°^W^x^°' 

where in the first integral x is taken on the first branch of x(£o) and in the second integral on the second branch. Let us make 
some natural assumptions. Initial distribution function is cut off at the typical energy of To, provided that To ^ 1 (bearing in 
mind that Tq as well as £ are in the units of energy £i (26)). Let us calculate N{£) in the region, where To ^ f -C 1, because 
it is seen from equations (36) and (36*) that in the case f ~ 1 the integral is small. Let us evaluate the contribution of the 
singularity x = 1/f in integrals (36) 

dx^ 
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Figure 1. The dependence of So versus x (the energy £ and the distance x are in the values of £i and r/Q (26) correspondently) 



Integral near the singularity is 



^^^^^sk^ J ^°(^°)T7#^r7?^'^^o + / /o(fo)4dfo. (37) 



1 f ^ v ^ 

■?2 / fo{£min) . = d£o ~ \/ £o ~ £min ■ 

^ J XOV to — tmin 

p 

Thus, its contribution is small. Because the main contribution to expression (30) is when x is far from we may use 

asymptotic values of x (see equation (27*)): x « \n{£ /£o) on the first branch; x l/£ + 1 on the second branch and obtain 
that 

£ oo 

'l/f-ln(fo/f)'^^°+f^ J 

£min £min 

Let us consider only £min To i.e. £ <C — l/lnTo because in the reverse case the integral is small. Now we can replace £min 
by and also 1/f S> ln(f /fo)- Equation (37) becomes 

s 

N{£) = ^ [ fo{£o){-ln£o)d£o + -r^, (38) 

tXo J t'^Xo 


where n = fo{£o) d£o. Rewriting expression (38) as 
we note that the first term is approximately equal to 

and is small compared to the second term because non-equality l/£^ S> — In To is the consequence of non-equality l/£ ^ 
— InTo. Therefore N{£) oc l/£^ and we obtain power law distribution with the index —3. The validity region of this formula 
To < f < — l/lnTo is quite large. 

We performed integration in equation (36) numerically having taken Gaussian initial distribution with different temper- 
atures. Results arc presented in Figs. 2 and 3. 

The value of power law index is minus three for small energies and rises slightly before the cut off of the distribution 
function. However, the value —3 of power law index is not applicable all the way down to f = 0. The natural boundary is the 
width of the acceleration region Ar (see expression (31*)). Because the main contribution to the integral (31*) is from the 
region near the point x = l/£ + 1 the lower limit for power law l/£^ is £i = 1/Ax. We also found mean power law index for 
the distribution of particles between energies £i and — 1/lnTo. Results are presented in Fig. 4. 

One remaining question in the formation of the spectrum of energetic particles is the origin of the initial particles, since in 
the first order on perturbed fields there is no acceleration of the cold plasma in equations (1). We have already mentioned that 
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Figure 2. Distribution function calculated with Gaussian initial distribution with temperature 2 • 10 (the energy S is in the values 

offi (26)) 



Figure 3. Distribution function calculated with Gaussian initial distribution with temperature 10 ^ (the energy £ is in the values of £i 
(26)) 



initial particles to be accelerated are originated in the region close to the Alfven resonance. It can be seen from equations (1) 
that in the second order on perturbed fields the acceleration always takes place and the power acquired corresponds to the term 
j • E in the energy conservation equation for the electromagnetic field. The damping of the wave occured due to the heating 
of the plasma by current ji which is described by second-order terms. If all energy of the wave is transmitted to particles, the 
condition /? > 1 would imply that the energy density of the wave is more than the energy density of the stationary field. Let 
this to be true at least on the surface of the resonance where the field is maximal. Then the energy acquired by one particle 
can be evaluated as To w Bq/Sttti. Given typical values Bo w 10~^G and n w O.lcm"^ we obtain To w lOMeV. 
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Figure 4. Averaged power law index (tlie initial particles energy To is in the values of £i (26)) 



4 SYNCHROTRON EMISSION. 

Let us now consider synchrotron emission of accelerated particles in a cylindrical jet with spiral magnetic field. Consistently 
with the results of section 3 we take isotropic distribution function of emitting particle in momentum and power law in energy 

dn = Ke£~''d£dVdnp (40) 

Here dn is the number of particles in the energy interval £,£ + d£, dV is elementary volume, dQp is the elementary solid angle 
in the direction of the particle momentum p, Ke — constant, 7 — constant. We will assume that ultralelativistic particles are 
distributed according (40) in the frame moving with plasma with the velocity 

u(r) = KB + Q'"(r)re^ 

We will consider two cases for the spatial distribution of particles, namely, uniform in the jet volume and concentrated close 

to the Alfven resonant surface r — ta- Influence of plasma on the omission of rolativistic c~ and e"*" is neglected. Absorption 
and reabsorption are also not taken into account. At first we assume that there is no Doppler boosting i.e. u = 0. In this case, 
according to e.g. Ginzburg, 1989, Stocks parameters are 

I=2±J^k{u) jiBsinxf dR, 

Q = k{i') j {BsinxY''^^'^^''^ cos2xdR, (41) 
U = k{v) I (Bsinx)'^^'''+'' 8m2xdR, 



V = 0, 

Here integration is to be fulfilled along the line of sight in the region filled with particles. 



3e 



1/2(7-1) ly- 

-1/2(7-1) £e 
D2' 



12 J mc^ 

X is the angle between the magnetic field B and the line of sight, x is the position angle in the picture plane of the electric 
field vector, D is the distance between the source and the observer. We will measure x clockwise from the direction parallel 
to the jet axis. For ultraxelativistic paxticles = 0, and the emission has only linear polarization. Degree of polarization of 
the outgoing radiation is expressed as H = \J Q'^ + U'^ / 1 , resultant position angle of the electric field is found from 

COS2Xres = , ^ , sin 2Xres = , ^ , < Xres < TT. 

Let h be the distance from the line of sight to the jet axis. We will assume that h is signed, i.e changing from —R to R for 
the different lines of sight. We will denote the angle between the jet axis and the direction to the observer as 5. In cylindrical 
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coordinates 



<P2 

7+1 ^ V sin6lsin20 ^ ^ 



n u ^ I \B±\'^'^''+'^ [Bl cos" ct> - {B. sing + B^sm4>cosef]h 

Q = k{^) / T^-Ti -0-2, d(i>, 43) 

' i>x smfc'sm <* 



7 IB, -h 
U = k{v) I ' ,^ 2S0Cos0(SzSin6i + B,^sin0cos0) — 



where = Trsgnft — arcsin (j>2 = arcsin It is seen that when one changes by tt — the expression for U changes its 

R, R. 

sign and, therefore, f/ = 0. Consequently, if Q > than Xres = 0, if Q < then Xres ~ 7r/2. Thus, the electric field vector 
either parallel or perpendicular to the jet axis. This fact may account for the bimodal distribution of observed Xr-es (Gabuzda 
et al., 1994a). In the case Bz = constant, B^ = Bzfl^r and 7 = 3 expressions for I and U from the unit jet length (i.e. (42) 
and (43) integrated over h) 

I = 2+lIlk{u) fdh [ \B^\^/W)f!__acl>, (44) 



7 + 1 



01 



Q = k{v) / dh / |Bx|V2(7-3) r^cos^0- (sing + n^rsin0cose)^]gg;^ ^^^^ 



-R <l>i 

may be reduced to the form 

R 

r „ ^9 „ 1 „ n 

dt, 



hl)/2^I_ [ , 



sine + Q," r{t){\ - -sin 9) 



n 

Q = fe(z/)B(^+i)/^7rsin6i j (a^V^(t) - 2)tdt. 





Prom above we see that in the case 7 = 3 the sign of Q does not depend on 9, but determined only by function Q,^{t). For 
the degree of polarization 

J{n^ r\t)-2)tdt 

TT 7 + 1 . 2 a 
n = -— sm 9- 



7 + 7/3 « 

/2idt(sin2 9 + n^V2(t)(l - i sin2 61)) 



In the case 7 = 3 and = 0(1 - (t/Rf) 

n = ^ii-J 3^in^^ 

4 



. 2 

sm 



,+^(,-l,i„'.) 



If |n| > 2\/6, then Q > 0, Xres = 0, i.e. orientation of the magnetic field in the jet derived from observations is perpendicular to 
the jet axis. If < |f2| < 2^/6 magnetic field looks parallel to the jet axis. Equations (44) and (45) were integrated numerically 
for = n(l — [t/R)^) and 7 = 2.1, which correspond to the power law index in the spectrum of synchrotron emission of 
(7 — l)/2 — 0.55, the mean observed value for cxtragalactic jets. Results are presented in Fig. 5. Positive values of the degree 
of polaxization correspond to Xres = 0, negative to Xres = tt/2. In the case 11 = radiation is entirely unpolaxized. Note also 
that in nonrelativistic consideration 11(180° — 9) = Ii{9) for any value of 9. 

In paper Istomin & Paxiev, 1996 we presented a hypothesis that knots observed in jets can be the manifestation of the 
"standing wave" phenomenon. If one knows the dispersion curves ui = u!{k) and wavenumbers femin at which du/dk = 
for perturbations one can calculate the velocity of moving crests of the "standing wave" as Vphas = i^min{kmin) / kmin and 
corresponding observed velocity Vvis = . Note also that changing sign of Ba, in (2) is equivalent to changing 

1 - Vphaa/C cos 9 
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Figure 5. Dependences of linear polarization 11 on angular rotational velocity Q (a) and on visual velocity of knots v^ig (b) for 
fundamental axisymmetric mode m = 0. 11 > when electric field vector is oriented parallel to the jet axis, 11 < when electric field 
vector is oriented perpendicular to the jet axis. Vyts > means that observed motion of knots is direct, v^is < means reversal motion. 
Crosses are for BL Lac objects and squares are for quasars listed in Table 1. 



sign of SI , since / and Q depend only on and and not on u or Er- On Fig. 5 we present the relationships of the 
polarization 11 of jet emission and v^is obtained in the case of il changing from to 100 for some angles 6. We also plotted 
observational points for jets in BL Lac objcts and quasars for which we were able to find both polarization and proper motion 
measurements in the literature and which are listed in Table 1. For each object we calculated degree of polarization and 
observed velocity as averaged over all bright knots observed in a particular object excluding a few knots peculiarity of which 
was indicated in original references. This pocedure might be very questionable although in the lack of observations of the 
emission from the space between knots on parsec scale we had nothing better. Numbers for Vvis were calculated assuming 
Hubble constant H — 100/ikm/s • Mpc 

In the case of accelerated particles concentrated very close to the Alfven surface r = va one can obtain from the 
expressions (44) and (45) the following formulae 



2'^ 



7 + 1 





" J sine 



6 + n^^(rA)f a(1 — sin'' 0sin^ ff) + 2D,'' rA{t) sin cos ^ sin 



7+1 



Q = k{v)B:,^ I sin^6l + n^^(rA)ri(l-sin^0sin2^) + 2n^(rA)r-Asin0cos6>sin^ 



^^O,^^ {rA)r A cos"^ <p — (sin ^ + (rA)rA sin ^ cos ^)^^ d<p. 



Since for axisymmetric modes (m = 0) always Iwminl > |fcmin| (Istomin & Pariev, 1994) one can see from the expression (5) 
for A that Alfven resonance does not exist for axisymmetric modes. Therefore, there should bo no particle acceleration by 
this mode. The next most important modes arc rn = 1 and m — —1. Since no one quantity kmin, ijJmin or 11 changes when m 
changes to —in and to — fi'^, it is enough to consider only m = 1 and all values of . Polarization 11 versus Vphas for this 
case is plotted in Fig. 6 as well as observational points. We showed first non-axisymmetric fundamental mode m = 1 which 
are expected to be excited more than higher modes with m > 1. Polarization curves are plotted only for values of at which 
Alfven resonance r = is present inside the jet. As previously, Q,^ = n(l — {r/Kf') and 7 — 2.1. It is seen that in the case 
of uniform particle distribution across the jet, curves for changing from 120° to 70° are located in the region where points 
from observations of BL Lac object (crosses) and quasars (squares) are situated. However, polarization curves for emitting 
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Table 1. Degree of polarization and velocities of knots in BL Lac objects and quasars 



Class 


Source 


n,% 




References 


BL Lac 


0454+844 


25 


0.6 


2,3 


BL Lac 


0735+178 


14.2 


5.5 


4,5 


BL Lac 


0954+658 


17.5 


5.7 


1 


BL Lac 


1803+784 


11.5 


1.8 


1 


BL Lac 


1823+568 


17.75 


4.05 


1 


BL Lac 


2007+777 


5.55 


2.9 


2,1 


quasar, 


1308+326 


3.77 


10.9 


6 


BL Lac? 










BL Lac 


0851+202 


6.2 


2.8 


7 


BL Lac 


1219+285 


30 


2.0 


7 


BL Lac 


2200+420 


10.9 


3.5 


7 


quasar 


0212+735 


-13.35 


3.9 


8 


quasar 


0836+710 


-31.33 


7.8 


8 


quasar 


3C279 


-15.37 


5.6 


8 


quasar 


3C345 


-7. 


7.7 


9 


quasar 


1928+738 


-9.15 


5.78 


8 


quasar 


3C454.3 


-15.55 


8.9 


8 



References.— (1) Gabuzda et al. (1994a), (2) Witzel et al. (1988), (3) Gabuzda ct al. (1989), (4) Baa th & Zhang (1991), (5) Gabuzda ct 
al. (1994b), (6) Gabuzda & Cawthorne (1993), (7) Gabuzda & Cawthorne (1996), (8) private communication of Gabuzda D.C., (9) Brown, 
Roberts, & Wardle (1994) 



Figure 6. The same plots as in Fig. 5 but emitting particles are concentrated close to the Alfven resonance surface, m = I fundamental 
mode. 



for the polarizatfon tensor Jag =< EaEg > and Stocks parameters Jag = \ [ jj ^ ^ ) when expressed in the source 



particles concentrated on the Alfven resonance surface, as it is the case if our acceleration mechanism works, do not match 
observation of BL Lac objects. Electric field vector is always oriented perpendicular to the jet axis. 

Now let us make calculations closer to the reality taking into account Doppler boosting and presence of strong electric 
field. We will denote values measured in the plasma rest frame as values with primes. Let us find the transformation formulae 

1 f I + Q U 
U 

frame XsVsZs and the observer frame XoUoZo- Using the transformation formulae for the electric and magnetic fields of the 
wave directly it is possible to obtain 

1 2/2 
(1 — 'u/c COS (5)^ 

Here we denote by S an angle between velocity u of the source and the direction towards observer k measured in the observer 
frame, correspondently S' is the same angle but measured in the source frame. Thus, polarization ellipse in the observer frame 
is obtained from polarization ellipse in the source frame by the rotation by the angle S — 5' in the plane containing the direction 
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to the observer n and u. Therefore the degree of polarization is not changed. Since for one particle Stocks parameters are 
expressed through the radiation flux with two main directions of polarization p^J} and p^} by the formulae / — + p^^^ , 

Q = — piy) cos2x , U = (pJV — p[^')sin2x then p^^ and p^^ are transformed as p^^' = p[3'^7:r— — 7—^— tt^i 

av (L — u/ccoao)-' 

c^f •* = p^'^} -j—^ — ^ . Therefore the intensity 

" dv [\ — u/c cos 5)2 

i{v,\^) = j {p^^){v',8',n,^,^')+p^^){v',£',n,^,,p'))x 

1~T-, 1 JvJ^et —--dtdQ-r'RdR, 

dv (1 — u/ccosoj^ dv 

where x' is an angle between momentum of the particle and the magnetic field vector, 1/)' is the angle between the direction of 
radiation and the cone formed by rotating velocity vector, both measured in the frame comoving with the fluid element. We 
will make integration over momenta in the plasma rest frame while integration over space coordinates in the observer frame. 
Integration over solid angle Qr' in the momenta space may be done analogous to the usual case of synchrotron emission of 
ultrar-relativistic particles, when there is no electric fleld E = — u x B. Integration over £' leads to (see, e.g., Ginzburg, 1989) 

i\ 7 + 7/3,, , f I 1 — u/ccosS\ . 1.2+1 dR 

/(.,k) = -^M^)y [^-^J^^j (B smx) ^ ^3^^ (47) 

1 — U / C COS S 

Here v' = f — — , S is the angle between u and the direction to the observer, v' is the angle between the magnetic 

field in the plasma rest frame B' and the line of sight, integration is performed along the line of sight. The expressions for Q 
and U are obtained in the same fashion. 

T+i cos2x dR 



^/ 1 \ 7 / \ /" / 1 — U/CCOSS \ ,r,l. l^2± 



TT/ ,\ ,/ \ f I 1 — m/ccos(5\ ^ ■ 1^2+1 sm2x dR 

U{v,k.) = k{v) / , ' (B'smx) ^ r 

7 1 Jl-uVc-^ / l-w/ccos( 



y^l -«2/c2 j l-u/ccosS' 
Using the transformation rules leads us to the expression for the B' sin x' 



sin^x' ^Bl + B'i-E'^- 



[—Er{l — u^)(uz sin ^ sin + u,^ cos 6) + u{cosd — u){u4,B^ + UzBz)^ 



(48) 



W*(l — W cos 5)2 

and using of (46) gives us the expression for x 

X = Xu + (T, (49) 
Er(v? — urns d){uz sin</)sinS + u^, cos 9) — s\v? 5{u^B^ + UzBz) 



tanx„ = 



ErV? sin 6 cos 0(1 — m cos 5) 

< X« < TT, 

ua, sin S cos 6 + Uz sin . ua, cos < 

cos a = : — ^ , sm (7 — — 



u sm u sm 

u cos 5 = — sin 9 sin (j)u^ + cos 9uz 
Assuming uniform particle distribution for the emission from the unit length we obtain 

, -1 + l/ik{v) ^ /l-wcosA""^ |B'sinx'l'^+'^/' 

/ = ^-^ rdr dff) ' 



7 + 1 sin 61 y VI - J 1 - M cos 5 



R r'Zn ■ - 



^ / J / ,,/l-ucos(5\ = , . /|(7+i)/2 cos2x 

"5 = -ti / / '^'^^ \B smx ^' z. ——T (50) 



sin 9 7y ' \^ ^/T^^n^ J 1 - w cos 5 

sm^X y„ '^\^yr^y' ' 1-ucos 
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Figure 7. Dependences of polarization 11 on Q (a), on for 7 = 2.1 (b) and for 7 = 3. (c). Emitting particles are concentrated close 
to the Alfven resonant surface. Here = +Q^rBz, fundamental eigenmode with m = 1. Notations are the same as in Fig. 5. In Fig. (a) 
solid curves are for 7 = 2.1, dashed are for 7 = 3. 



Assuming particles to be located close to the Alfven resonance surface in the region |r — r^l < A and A <C i? we have 



^ 7 + 7/3fc(j/) ^ fl--ucosS\ ^ |B'sinxf^+'^/^ 

'=TTrih^^/ '^^'^-^T^) l-.cosJ 



^ ^ A / d<t>rA ' 



sine Jo \ VI -M^ 

|B'sinxf"+^>/^ g°«2x (5^) 
1 — It cos d 

fc(j^) . f^^ ,, /l— «cos(5 



sine 7g V ^ 

'sinx 



D/ • 'i(7+i)/2 sin2x 
B smx : 



1 — M cos 6 ' 

In the force-free approach constant K determining the component of u parallel to the magnetic field in the observer frame 
B remains to be free. We choose it such that the velocity u of plasma has the least possible absolute value for a given Bz 

B B B^ 

and Q^r. This choice is the same as in (3) Uz = — „ — ^fi^r, — —-^ — ^^^^^ value of K equals to the 

Bfi) ^" Bz B^ -\~ Bz 

—Q^rB^/B'^ . It can be checked that under change of (p to —cf) the value of Q is not changed, and the sign of U is reversed. 
Thus, as a result of integration U = and either Xres = or Xres = 7r/2 as it is when we do not take into account Doppler 
boosting and electric field. 

We performed numerical calculations according to (50) and (51) using (48) and (49). As above we take Bz = constant, 
= f2(l — {r/Rf), B^ = ±Q^rBz. Results for uniform particle distribution do not match observations. Taking into account 
presence of the electric field and Doppler boosting changes the polarization of synchrotron emission essentially. Now only 
for small values of angle 6 and high Q apparent direction of the magnetic field as it is usually derived from polarization 
measurements (just perpendicular to the electric vector in the polarized component of the radiation) is transversal to the jet 
axis. For most of inclination angles 9 and angular velocities f2 apparent magnetic field is longitudinal. One can also see from 
Fig. 7, where we plotted 11 versus Vvis for B^ = +Q,^rBz, that our simulation does not match observational data. However, 
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Figure 8. Dependences of polarization 11 on Q (a), on v^is for 7 = 2.1 (b) and for 7 = 3. (c). Emitting particles arc concentrated close 
to the Alfven resonant surface. -B^ = —Q^rBz, fundamental eigenmode with m = 1 Notations are the same as in Fig. 5. In Fig. (a) solid 
curves are for 7 = 2.1, dashed are for 7 = 3. 



for B4, = -D.'^rB, wo arc able to reproduce observations well enough, if particles are distributed close to the Alfven surface 
and m = lorm = — 1. As for the case without Doppler boosting ofTcct and without taking account electric field, n is not 
changed when one changes the sign of Q,, but 11 is changed when one put = —Q,^rBz instead of B^ = Q,^rBz. Let us list 
all transformations which do not change equations for radial modes (5) or lead to complex conjugated equations for 

i) uj ^ —uj*. k —k, Q —Q; 

ii) Lu — > —Lu*, k —^ —k, m —> — m; 

iii) n — » — f2, m — » — m; 

iv) k -k, B^ = +n''rBz ^ B^ = -Q^rB^; 

v) Bz -Bz. 

So, if one set of parameters satisfies equations (5) than parameters with changed signs will also satisfy equations (5) for the 
same ^ or for ^* . Obviously, these substitutions are valid for kmin and ujmin as well. The value of rA also is not changed under 
transformations listed above. Degree of polarization D is not influenced by changing sign of Bz as well as by changing sign 
of n. The crucial fact is that changing sign in relation = ±Q.^rBz leads to changing sign of Vphas = oJmin/kmin, whilst all 
other four transformations change neither 11 nor they change Vphas- Therefore, in the sense of correlation between 11 and Vvis, 
which can be verified experimentally, all different cases are reduced to considering both signs in the relation B4, — ±fl^rBz, 
but we can restrict ourself to Bz > 0, m > and all O for which va is located inside the jet. Using formulae (48), (49), (50) 
and expressions (3) for Uz and one can check that Stocks parameters remain unchanged when one changes sign of B^ and 
replace angle 6 hy n — 6. Also, simultaneous changing of sign of tiphas and 9 to n — 6 leads to that Vvis is transformed into 
—Vvis- Therefore, dependences of 11 on Vvis in the case can be obtained from same dependences in the case 

B^ = +Q,^rBz by reflection Vvis —Vvis and relabelling each curve from 9 to n — 9. 

We plotted degree of linear polarization fl versus Vyis for m = 1 in Figs. 7 and 8 for B^, = +n^rBz and B4, = -fl'^rBz 
correspondently. Alfven resonance exists for 0.7 < Q < 4.9. Therefore, along curves in Figs. 7 and 8 the value of Q changes 
in this interval. We see that curves in Fig. 7 do not match observational data for knots in jets, while those in Fig. 8 do. 
Dependence on the power law index 7 of emitting particles is small. We plotted results for 7 = 2.1 and 7 = 3. which are 
maximum and minimum values of 7 obtained in the process of the formation of the spectrum (see section 3). 

We can suggest the following explanation of the fact that knots in jets with B^ = —OFrBz are observed more often than 
in jets having the opposite winding of the magnetic field. If we look at the models of the formation of magnetically driven jets 
from accretion discs (see, e.g. Pelletier & Pudritz, 1992) we will find that magnetic field lines are always bent in the direction 
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opposite to the direction of disc rotation as soon as the velocity of the outflow is directed outwards from the disc plane. This 

means that in the jet approaching the observer = —fl^rBz, whilst in the jet receding from the observer = + fl^rBz. 
Disturbances propagating along approaching jet should be brighter than those in the counter-jet due to the Doppler boosting 
effect. Therefore, we will observe knots in parts of jets having = —^^rBz more often than in those with B^ = +Q,^rBz. 
This explains the fact that the number of sources having 11 > is approximately the same as the number of sources having 
n < despite that for most of the values of parameters calculations give longitudinal polarization 11 < 0. Also the difference 
between jets in BL Lac objects and jets in quasars finds its explanation by systematically different values of . In the frame 
of our model jets in BL Lac objects have higher value of Q.^ than jets in quasars. Although we performed calculations only for 
one particular profile of which seems to be close to what one should obtain from theory of the formation of jets (Istomin 

& Pariev, 1996) we found this tendency and it should be general for profiles of which are qualitively similar to chosen 

one. There are many physical arguments though (Istomin & Pariev, 1994; Istomin & Pariev, 1996; Beskin, Istomin & Pariev, 
1992a), why the real dependence f2^(r) should have qualitative properties similar to 0.(1 — {r/Rf) that we chose in our 
computations. 



5 SUMMARY 

In present work we considered possible acceleration of electrons and positrons inside relativistic rotating electron-positron 
force-free cylindrical jet with spiral magnetic field. It is very plausible that inner parsec-scale jets in active galaxies can be 
described in the frame of this physical model. The observations of jots in 3C273 and M87 points out that the cxtragalactic 
jets are most probably electron-positron rather than electron-proton ( Morrison & Sadun, 1992; Reynolds et al., 1996 ) Also it 
seems to be very likely that such jets are not in steady state and contain waves excited by environmental processes, of which 
most powerful would be the variability of accretion rate onto central black hole and accretion disc instabilities both strongly 
perturbing magnetic field at the base of the jet. We considered the behaviour of such excitations inside cylindrical force-free 
jet in our previous works (Istomin & Pariev, 1994; Istomin & Pariev, 1996) and found that in wide range of parameters 
determining the equilibrium structure of electromagnetic fields and perturbations inside the jet there exist resonant surfaces 
on which the phase speed of eigenmodes is equal to local Alfven velocity. Those eigenmodes for which there is an Alfven 
resonance occured to have small damping decrement, whilst those having no Alfven resonances are neutrally stable. This fact 
points that existence of Alfven resonance leads to losses of energy of corresponding eigenmodes. In present work we found how 
particles are accelerated in the strong fields near the Alfven surface. Acceleration process and synchrotron losses combined 
together form power law energy spectrum of ultrarrelativistic electrons and positrons with index between 2 and 3 depending 
upon initial energy of injected particles. The power law spectrum goes up to the energy Emax, where the sharp fall down 
occurs (see Figs. 2,3). The magnitude of Emax depends on the initial particle temperature To as well as on the characteristic 
acceleration energy fi(26), £max ^ fi/ln(fi/To). The quantity of Tb evaluated from the equipartition condition in the 
acceleration region is equal to 10 MeV for the quantities of magnetic field and density 10~'^G and O.lcm^"' correspondingly. 
As to the quantity of £i, it is, in accordance with (26), of the order of 1Q^{5B / BofUiN for the same value of Bo, where SB /Bo 
is the dimensionless amplitude of the perturbation. For large perturbations [SB w So) £i reaches the values of the order of 
lO^MeV. Thus particles accelerated near Alfven resonance are in the energy range of lOMeV < £ < lO^MeV. These particles 
emitt synchrotron radiation in the range of frequences approximately from IMHz to lOOGHz which covers the frequences of 
modern radio observations. 

In order to obtain observable predictions of our theory of standing wave excitations inside jets, we computed synchrotron 
radiation emitted by accelerated particles. We considered two possibilities: 1) ultra-relativistic particles distributed uniformly 
across the jet, 2) ultra-relativistic particles concentrated close to the Alfven resonance surface. The latter distribution follows 
from our acceleration theory and wo compared it with the former one which might be produced by other acceleration processes. 
We assumed uniform magnetic field in the jet, so the equilibrium configuration of the jet is determined only by the radial 
profile of angular rotational velocity f2^(r') of magnetic field lines. In our computations we also specified component of 
the velocity parallel to the magnetic field lines in a way to minimize the kinetic energy associated with the flow and took 
r2^(r) = n(l — (r/Rjet)^)- The actual velocity is determined by the processes of jet formation which is poorly understood 
now, still such minimal velocity is physically favoured (Novikov & Frolov, 1986). The prominent characteristic of synchrotron 
radiation is its polarization: degree of linear polarization and orientation of electric vector with respect to the jet axis. We 
calculated these quantities for the emission integrated over the jet cross section taking into account relativistic bulk spiral 
flow of plasma in the jet, which leads to the presence of electric field comparable by magnitude with the magnetic field, 
Doppler boosting effect for the intensity of synchrotron radiation and swing of the polarization position angle. We want to 
stress that the last effect, which was first pointed out for the simpler case of uniform magnetic field parallel to the jet axis 
by Blandford & Konigl in 1979, substantially modifies resulting degree of linear polarization from the whole jet. However, 
because of the axial symmetry of the problem, only two orientations of polarization position angle are possible: parallel or 
perpendicular with respect to the jet axis. Observations of polarization of VLBI scale jets (Gabuzda et al., 1994a) also show 
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bimodal distribution of the position angle: most of knots in BL Lac objects and quasars have polarization either close to the 

jet axis or perpendicular to it. 

We also calculated velocities of crests of standing eigenmodes and compared them with the observed proper motions of 
bright VLBI knots in jets. We were able to match observational data only for the case when emitting particles are concentrated 
close to the Alfvcn resonance surface and magnetic field lines are twisted in the direction opposite to the direction of rotation 
n^. The last fact can be understood by selection effect due to the Doppler boosting which makes sources with jets approaching 
the observer more favourable observed than receding jets. Comparing calculations with observations we can estimate angular 
rotational velocity of jets. We obtained that Q.^ in BL Lae objects is intrinsically larger than in quasars which accounts for 
the difference in polarization properties of these two classes of objects (Gabuzda et al., 1994a): in BL Lac the magnetic field 
inferred from polarization position angle is oriented predominantly perpendicular to the jet axis, in quasars it is parallel to 
the jet axis. 
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